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We discuss how new data for = (g — 2)^/2 improve constraints on new physics. Using two 
types of estimations of ajj"'' (Davier & Hooker (case A) and Jegerlehner2000 (case B)) we evaluate 
95% CL intervals for a new contribution, which can be used to constrain any model beyond the 
Standard Model. We apply these intervals to the general 2HDM ("Model 11"), where up to now one 
light neutral Higgs particle, h or A, was allowed by the data. Assuming that only one Higgs boson 
contributes to af/^^'^' the two-loop calculation based on the case A leads to the exclusion of scalar 
h while a pseudoscalar with mass between 10 and 70 GeV ( for tan/3 > 20 ) is allowed. For case B 
the upper limits for tan (3 for both scalar and pseudoscalar are obtained. 
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I. INTRODUCTION 

The precision measurement of g — 2 for the muon is 
expected to shed light on "new physics". The new E821 
result based on 1999 data [|| leads to a current mean of 
experimental results for {g — 2)^ (from 0|) 

^exp ^ ^ V ^ (15) . 10-10, 

where the accuracy of this result (in parentheses) ap- 
proaches the size of electro- weak contribution, a^^. The 
ultimate accuracy of the E821 experiment is 4 ■10^^'^. 

The QED and EW contributions to a^*^ are well under 
control. The predictions for the hadronic contribution 
ijhad^ which is more than forty times larger than a^^, 
differ considerably among themselves both for the central 
value and its uncertainty. This uncertainty is presently 
of order (7— 10) • 10"!", with the dominant error coming 
from the vacuum polarization contribution. Useful 
discussion of various estimations of a'^'^ can be found in 
lH , see also more recent papers 0,^, ^, @, ^, |§1 . 

The difference between the experimental data, a^fP, 
and the Standard Model (SM) prediction [2-47], a^*^, 
defines the room for " new physics" . Obviously the un- 
certainties of the hadronic contributions influence the es- 
timation of a size of new effects. To illustrate the present 
situation we calculate 95% CL intervals {lim{95%)) for 
an allowed new contribution, Sa^^, using two representa- 
tive SM predictions : one based on calculation of a'^"''^ 
with the vacuum polarization estimation by Davier 
and Hocker |2^, Q and the other by Jegerlehner ^ . 
We derive also intervals which may be relevant for models 
leading to a positive contribution only (L'm+(95%)). 

The obtained intervals we apply to constrain the 
parameters of the CP conserving 2HDM ('Model IF) 
p8| , This model is based on the two doublets of 



complex scalar fields, and predicts existence of five Higgs 
particles: two neutral Higgs scalars h and H, one neu- 
tral pseudoscalar A, and a pair of charged Higgses H^. 
In MSSM, which has a Model II Higgs sector, the mass 
limits (95%CL) from LEP experiments are, e.g. for 
a maximal Mh scan: M/j(M^) greater than 90.0 (91.9) 
GeV and tan/3 < 0.5 or tan/3 > 2.4 l^]. The present 95 
% CL hmit on the SM Higgs mass is: Mh > 114.1 GeV 
!§. The SM Higgs particle with mass around 115 GeV 
(corresponding to a maximum likelihood ratio from LEP 
data 1521) contributes to the a^*^ at the level of 10-^^ 
(two-loop result ||45[|). 

In the non-supersymmetric 2HDM (II), which we study 
here, one light neutral Higgs boson h or A with mass be- 
low 50 GeV is still allowed by data [50-77], see e.g. the 
LEP results from Higgs-strahlung and hA pair produc- 
tion [|5j, [52), from the Yukawa processes Z — > ffh{A) 
with / = T or 6 |^ , p5[ and the process Z h(A)j 
[|§, see also |5^, The dedicated fit to the EW pre- 
cision LEP data performed within this model, even with 
masses of ft. or A below 20 GeV, is equally good as the 
corresponding fit in the SM or MSSM [|9). Unfortu- 
nately a potential of the HERA collider, discussed in 
[]60| , |6l| , has never been explored to put limits on very 
light Higgs bosons of the 2HDM (II) . Future Linear Col- 
lider experiments are not expected either to close a light 
Higgs window in the 2HDM (II) [^ , even these planned 
to run with a very high luminosity, once more at the Z- 
peak (GigaZ) |6^. On the other hand one should keep in 
mind that some theoretical arguments disfavor the Higgs 
scalar h with mass below 90 GeV in 2HDM(II) Q. 

A light Higgs 'could conceivable evade discovery at 
LEP and yet show up in a analysis of a low energy data' 
as {g — 2)^, as it was pointed out in paper [Q. We 
have used previous experimental data and the SM pre- 
diction(s) to constrain 2HDM (II) in [Q. That one-loop 
analysis led to a very small improvement in comparison 
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to LEP limits. It was expected that with increasing pre- 
cision the (g — 2)^ measurements would lead to more 
stringent constraints. As we have pointed out in [ |66| , 
the sign of the one-loop 2HDM (II) contributions to 
is correlated with the type of the lightest particle: h, A 
or H^. We have described the condition which would 
lead to the exclusion at 95 % CL of a light h or A. This 
condition is presently fulfilled for a pseudoscalar if the 
DH result for aji"'' is used in the one-loop approach, see 
also@ However, as it was pointed out in 
the two- loop calculation leads to very different results. 

In this study we constrain 2HDM (II) by the new data 
[01 using the two-loop calculation. We apply two dif- 
ferent hadronic contributions a';"'' as mentioned above 
(DII=case A and J2000=case B) and two types of 95% 
CL intervals (lim and lim,-^). Obtained constraints we 
combine with constraints from other processes. 

The one- loop analysis of the 2IIDM(II) with results 
which partly overlap with results of this paper can be 
found in |Q, see also The two- loop analysis, which 
we follow here, is given in |8|] . Related study can 
be found in ||82|] . Studies of g — 2 for muon within the 
context of supersymmetric models were performed ear- 
lier js^ , many new analyses have appeared recently 
Dark matter problem in supersymmetric 
Relations to neutrino masses 




models is discussed in [89 



and mixing are studied within supersymmetric and non- 
supersymmetric models in |Q . Other new analyses made 
recently are [^, analyses based on the general 2IIDM 
(Model III) are presented in Analyses of some 

aspects of cosmic rays can be found in Relations 
between the electric and magnetic dipole moments are 
studied in |9^ . An model-independent analysis was done 
in (9|. 

In Sec. II we discuss the new {g — 2)^ data and de- 
rive limits for a new contribution which can be used in 
any model beyond the SM. In Sec. Ill we apply the ob- 
tained limits to 2IIDM (II) and derive constraints on the 
parameters of the model. Next we combine these con- 
straints with other experimental information. We study 
separately the one-loop and two-loop results. Sec. IV 
contains conclusions. 



The QED results can be found in |]T^, |9j . Hadronic con- 
tributions were obtained in: the leading vacuum polar- 
i^ation tcrmfv^l) |, |J,Jo[ 12^, ||, g ||, §[ |§, 
P^ , pg] , |30| , |3l| , |32[ and ^ST^J the higher order vacuum 
polarization term (vj3.2) |2C[ P2| . [2q , pSj, t he light-on- 
light term (Ibl) §, 
results are given in 



12^, 13 




J. The EW 
J||,E|,|7). 

A useful compilation of recent results is given in see 
also 



00- 



The error of the hadronic contribution dominates the 
total error of the SM predictions. Moreover the hadronic 
contributions calculated by various authors as discussed 
in ||, ^, ||, ^, 1^, ^, ^ differ strongly leading to distinct 
estimations of a size of the SM contribution and therefore 
also of the 'beyond SM' effect. All of these estimations 
can be divided into two classes, depending on whether the 
SM prediction for {g — 2)^ is or is not in agreement with 
the data. We think that it is sensible to study separately 
consequences of these two classes of the SM predictions. 
It allows to illustrate the present sensitivity to ajj'"', (see 
|66| for previous results based on a similar approach,) 
moreover it may serve as a guide for the future results, 
which obviously will belong to one of these two classes. 

In the analysis we consider two representative cases: 
case A based on Davier and Hocker calculation of the 
leading vacuum polarization diagram (v.p.l) [ p9| , pO[ and 
case B based on the corresponding Jegerlehner calcula- 
tion j26[ ^ , with a smaller and larger hadronic contribu- 
tion (and its uncertainty), respectively. In the DH analy- 
sis the e^e~ and r-decay data are used, while J2000 uses 
only the e^e~ data. The higher order hadronic (v. p. 2) 
contribution and the light-on-light contribution are taken 
from [22, 3^ and [37[ 3g, respectively. In the table below 
we collect, following |^the corresponding SM contribu- 
tions (and their uncertainties). 



case 



A [in 10" 



B [in 10 



-111 



QED 116 584 706 (3) 
had 6 739 (67) 
EW 152 (4) 



116 584 706 (3) 
6 803 (114) 
152 (4) 



tot 



116 591 597 (67) 116 591 661 (114) 



II. NEW G-2 DATA FOR THE MUON 

A. New [g — 2)^ results 

The current world average experimental data on [g — 
2) for muon averaged over the sign of the muon electric 
charge is given by ( from [Q): 

a^i^^ = 116 592 023 • 10"", with cr^xp = 151 • 10"". 

The Standard Model prediction for this quantity con- 
sists of the QED, hadronic and EW contributions: 



The difference between two predictions of the SM, case 
A and case B, is 64 • 10~^^, with the corresponding differ- 
ence in accuracies equal to 47 • 10~^^. It is worthwhile to 
compare these numbers to the ultimate accuracy of the 
E821 experiment 40 • 10"". 

Results based on other estimations of the hadronic con- 
tribution, can be easily obtained from results for consid- 
ered cases, see discussion below. 

B. The room for new physics 

Here we present results of calculation of the difference 
between the experimental and theoretical SM results for 
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(g — 2)^, which can be used for any model going beyond 
the SM. First we derive Aa^, equal to the difference of 



the central values, a'^^ 



a„*^ = Aa^, and the error a 



for this quantity. Knowing Aa^ and a one can calculate 
an allowed, at chosen confidence level (CL), interval of 
an additional contribution. For the purpose of this study 
it is enough to obtain a by adding in quadrature the 
corresponding experimental and theoretical errors {aexp 
and atot) and to assume a Gaussian distribution. Under 
these assumptions we calculate in both cases, A and B, 
the 5a regions, symmetric around Aa^, allowed at 95% 
CL. 



case 


A [in 10-"] 


B [in 10""] 


Aa^(cr) 


426(165) 


362(189) 


lim{95%) 


102 < Sa^ < 750 


- 8.65 < Suf, < 733 


positive Sttfj, 


99.5% 


97.2% 


lim+{95%) 


109 < Sttf, < 744 


28.5 < Sttf, < 696 



In the standard approach one obtains intervals called 
lim{95%), see the above table for results. We see that 
although at the one sigma level, i.e. for the interval 
Aa^ ± a, the allowed additional contribution to a^ in 
cases A and B are of a positive sign only, at the 2 a 
level or 95 %CL the more conservative estimation of the 
hadronic contribution to a^^^ (case B) makes the nega- 
tive possible. The SM prediction lies within the 95% 
CL interval for case B, while for case A it is outside the 
corresponding interval. 

This difference leads to different forms of the limits for 
a new contribution. The 95 %CL interval leads in case 
A to an allowed positive contribution (an allowed band) 
and at the same time to the exclusion of the negative con- 
tribution (at the higher CL level, see below). For the case 
B, the positive (negative) contribution is only bounded 
from above (below) {upper limits for the absolute value 
of the new contribution). That means that presently the 
accuracy of the theoretical predictions for the hadronic 
contribution (both ajj'"' and ahad matter !) influences in 
qualitative way the constraints on the new physics ||9g[| . 

Now we discuss consequences of the present (g — 2)^ 
data for models which can give contribution of only one 
sign. We see (the table above) that for both A and B 
cases a negative Ja^ contribution is very unlikely: a pos- 
itive (negative) contribution corresponds to 99.5 (0.5) % 
CL for A, while for case B to 97.2 (2.8) % CL. At this 
level models leading to only negative Sa^ can be excluded 
or saying differently, within models which give a definite 
sign contribution to only these which give a positive 



6a^ can be realized at 95% (or higher) CL. 

To obtain the allowed range of parameters of such a 
model we calculate the corresponding 95% CL intervals 
normalized to the positive contributions only |97|. Re- 
sults for these intervals, called lim+(^%), are presented 
in the last row of the table above [100|. Of course, these 
two types of 95 % CL intervals. Urn or lim+, will lead 
to very similar constraints of parameters of the model 
for case A, since they will be obtained from only slightly 



different, due to the 0.5% change in normalization, al- 
lowed 6a^. However, significantly different constraints 
will arise for case B. Instead of upper limits obtained 
in the L'm(95%) approach, in the L'm+(95%) method an 
allowed band for a positive contribution is obtained. 

We observe that the maximal positive 6a^ values dif- 
fer less than 10 % for all discussed cases (696 to 750 
in 10""'^^). The minimal positive i5a^ obtained using 
lim+{95%) method differ by a factor 3 in cases A and 
B, 109 and 28.5 in 10"", respectively. 

Other recent estimations of the hadronic contributions 
lead to the following Aap((T) and intervals (5a^ in [10^^^]: 



Reference 



Hm(95%) 



Jegerlehner[6] J2001 376(186) 12 < Sa^, < 740 
Narison [5] N 375(170) 41.8 < Sa^ < 708 
Melnikov [7] M 377(216) -47.2 < < 801 

DeTroconiz 

-Yndurain[8] TYl 363(184) 2.52 < Ja^ < 723 

TY2 338(171) 3.08 < Sa^ < 673 

Prades [9] P 403(169) 71.8 < Sa^ < 734 

Note similarities in the obtained Aa^ values and at the 
same time large differences in e stimated uncertainties a 
for the first three analyses |101|. The first two analyses 
are similar to case A, while the third one has properties 
of the case B. The maximal Sa,, differ from the corre- 



sponding numbers for cases A and B within 10% , the 
minimal (5a^ differ much more. 

The results of the newer analyses also are presented in 
the above table. The TYl analysis which corresponds to 
a (v.p.l) calculation based on the (new) e^e~ data, and 
TY2 where both the e~^e~ and r decay data are included 
lH, the P results are based on the weighted average (for 
v.p.l) of the averaged estimation of (J2001 and TYl) 
and (DH and TY2) Note that all these estimations 
of a'^"''^ are based on a Chiral Model for the light-by- 
light contribution |3^ . The Quark Model gives very 
different results. The corresponding limits were derived 
in g: the estimations TY3 and TY4 (analogous to TYl 
and TY2, i.e. without and with r decay data) leading 
to -177 < 6af, < 547 and -177< fc^ < 497 in [10""], 
respectively. 

In the following analysis we will apply intervals ob- 
tained for case A and B, a simple rescaling allows to 
translate the final constraints to results relevant for any 
other present or future estimations of i5a^ . 



III. CONSTRAINING THE 2HDM (II) 

A. A model 

In the non-supersymmetric CP conserving 2HDM the 
Higgs sector contains the two neutral scalars, h and H, 
pseudoscalar A and charged Higgs bosons H^. Beside 
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their masses, three parameters: tan /5, which is the ratio 
of the vacuum expectation vahies of the Higgs doublets 
W2/W1, the mixing angle in the neutral Higgs sector a, 
and in addition one more parameter, e.g. the Higgs 
self-coupling ghH+H- 1 specify the model. 

In the Model (II) implementation of the 2HDM, one 
doublet of fundamental scalar fields couples to the u- 
type quarks, the other to the d-type quarks and charged 
leptons (this way FCNC processes are avoided at the tree 
level) [Q. The ratios, relative to the SM values, of the 
direct coupling constants of the Higgs boson h or H to 
the massive gauge bosons V — W or Z , and to fermions 
( i.e. Yukawa couplings) can be determined via angles a 
and /3 . For Xi = 9il{9'l)sM (and similarly for H) 

we have, in form suitable for discussion simultaneously of 
h and H: 



sin(/3 — a) , = cos(/3 — a), 

=X^+cot/3xy, = Xy -cot/3x^ 

Xd =X^-tan/3xv, xf = Xv + tan/3x^, 



(1) 
(2) 
(3) 



with (xy)^ + (Xy)^ = A very useful pattern rela- 

tion among these couplings holds for both h and H M : 



iXu + Xd)xv = 1 + XuXd- 



(4) 



For Xv — 1 couplings of h have the SM values, cou- 
plings of H to gauge bosons are equal to zero while cou- 
plings of H to fermions may differ co nsid er ably from 
the SM values, for smaU or large tan/3 |l0|, (lo|. For 



Xv = 1 the H is SM-like while h has different properties 
(1-3), e.g. Xd be very large for large tan/3. From 
(2) and (3) or (4) it follows that for xv = one obtains 
XuXd = -1- 

For the pseudoscalar there is no coupling to W/Z. 
The Yukawa couplings to fermions XdiXu) contain tan/3 
(cot /3) factor. The Xd large (small) for large (small) 
tan/3 value, respectively, with the opposite pattern for 

Xu ' 

In this analysis the Yukawa coupling Xd, relevant for 
a Higgs boson coupling to a muon, plays a basic role. 
It is equal to tan /3 for a pseudoscalar and and, if 
in addition xv — sin(/3 — a)=:0, also for a scalar (more 
precisely Xd = ±tan/3). 

At the two-loop level there appears a possibility of hav- 
ing a charged Higgs boson in a loop (see below). A cou- 
pling of _ff + to a scalar h has a form: 



x'h. 



(1 



Ml 



iA+x'u). (5) 



with the normalization as for an elementary charged 
scalar particle in the SM. For Xv ~ "^'^^ S^^s 



Xli. 



M, 



2 _ ,,2 A,f2 _ 2 1 



B. Existing constraints 

Main constraints of the parameters of the 2HDM (II ) 
come from LEP experiments, see also a discussion in [ |67[ . 
From the Z — > Zh process upper limits on sin^(/3 — a) 
were derived ||, ||. Fr om the tightest limits it 
follows that sin (/? — a) should be smaller than 0.1 for 

the O^Mh^ 50 GeV, and even below 0.01 for a lighter 
scalar. The data for the cross section for the pair (ft,. A) 
production, proportional to cos^(/3 — a), when combined 
with the Z Zh data leads to an exclusion of a small 
mass region m the {Mh,MA) plane[|l|, ||. According 
to these data the 2HDM (II) may accommodate a very 
light ( ^ 45 GeV) neutral scalar h, with small sin(/3 — 
a), or a very light pseudoscalar A as long as: Mh + 
> 



50 |51 



Ma 

A neutral 



(6) 



5^, see also the newest results ]58| ]. 
Higgs particle at LEP I has also been 
searched for in the Yukawa process, e+e~ — > ffh{A), 
where / means here b quark or r lepton. For a light 
scalar this is an additional, and if /? = a, the most impor- 
tant source of information. A (still preliminary) ALEPH 
analysis of the Yukawa process for a pseudoscalar |Q led 
to the exclusion at 95% CL for the (tan /3, Ma) plane, al- 
lowing for a large tan/3, above 20, for mass larger than 
2 GeV. Similar analysis was performed by DELPHI for 
the b quark-couplings to a scalar h and pseudoscalar A. 
New measurements of Yukawa process at LEP by OPAL 
and DELPHI groups |5^, improve slightly these results, 
see discussion in Sec. IV. 

A measurements oi Z ^ h{A)j performed by all ex- 
perimental groups at LEP I was used to obtain upper 
limits (however weaker than from the Yukawa processes) 
and lower limits on the Yukawa couplings Xd p6| . Still 
large part of the parameter space remains unconstrained. 

Also the dedicated global fit to the EW precision data 
in the 2HDM (II) framework allows for an existence of 
very light scalar or pseudoscalar | ]59| , for a partly con- 
strained the heavy Higgs boson sector, including the H 
and bosons. Note that the lower mass limit of 
estimated from the direct search at LEP is 78.6 GeV [ pO[ . 
The 6 — ^ 57 data interpreted in the 2HDM (II) give, ac- 
cording to the newest results ||6^, M/j± > 320 GeV or 
even higher ||69|| . One of the important message from the 
global fit based on the EW precision data |5^ is that for 
a light h and large tan/3, and for mass of Mh below 1 
TeV, an upper limit for Mjj± can be derived. In order to 
agree with the above lower mass limit based on b —^ sj 
data, tan /3 should be smaller than 22 (28) for a light h 
with mass 10(20) GeV. For the upper mass of H equal 
to 5 TeV, these maximal value of tan/3 increases by ~ 3. 

In light of the above results one can conclude that there 
is still a possibility of the existence of one light neutral 
Higgs particle h or A with mass even below ~ 40-50 
GeV. Since for a very light h the limit sin^(/3 — a) ^ 1 
should hold, the second scalar, H, is expected to mimic 
the SM Higgs boson couplings, as discussed in Sec. Ill A. 
Therefore it is reasonable to expect that its mass is equal 



5 



to '--^ 115 GeV (or slightly higher). 



Other low energy measurements do not contradict such 
scenarios. The 77 decay data exclude only a very 
light h, with mas s be low 280 MeV. The Wilczek process, 
T — !■ h{ A)j I74 |104|, has been measured by few groups 
Jt^ , [75 , 76| 77|]. Unfortunately the corresponding predic- 
tions have large theoretical uncertainties both due to the 
QCD and relativistic corrections, see [Esl 73 . In addi- 



tion, in some experimental analyses the production of the 
Higgs boson was treated according to the general 2HDM 
while for a decay of such Higgs boson the SM rates were 
assumed, e.g. |7^. All of these measurements of the T 
decays give only upper limits for a coupling of the Higgs 
boson to b quark. Although for mass below 10 GeV these 
limits seem to be stronger than others mentioned above, 
large uncertainties of various sources make difficult to use 
these limits on a similar footing as the LEP ones. In the 
present analysis on 5 — 2 data , in which as we will see 
below also lower limits for the Yukawa coupling appear, 
the data for the Wilczek process, even with large uncer- 
tainties, will play an important role in closing low mass 
(below 10 GeV) part of a parameter space for 2HDM (II) . 
In this analysis we apply three different constraints from 
the T h{A)j process, from (denoted in figures as 
K), ^ (N) and (L). 
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Unfortunately there are no limits from the HERA col- 
lider on very light Higgs bosons [^0|, |6l|. There are im- 
portant upper limits for the Yukawa couplings Xd for h 
and A from the TEVATRON data in the large mass re- 
gion. They were obtained originally for the Higgs bosons 
in MSSM M: for mass say 70 (120) GeV tan/3 should 



be above 34(60) |70 . These limits, rescaled by a factor 
•\/2, should hold also for h and A in context of the 2HDM 
(II). 

In Fig. 1 we present the existing 95% GL limits for 
a Yukawa coupling Xd for h (solid line) and A (dashed 
line) in the 2HDM (II). Upper limits are from the Yukawa 
process (data from ALEPH ||| for A and DELPHI ^ 
for h and A). Lower limits from the Z h{A)'^ processes 
p6| measured at LEP are also shown. In addition the 
upper 90% CL constraints from the T decay from (K) 
are shown, rescaled by a factor 2 to take into account the 
preference of the Higgs boson to decay into tau leptons 
(It occurs in 2HDM (II) for the considered mass range, 
for Xd bigger than 2.). The TEVATRON limits [|o) (with 
a factor y^) are displayed as well. 

The one-loop analysis based on the previous .g — 2 for 
the muon results |66| led to a slight improvement of the 
upper limits of Xd as compared to the Yukawa processes, 
for Mfi(A) < 2 GeV (not shown). The new data improve 
constraints of the 2HDM(II) considerably, what will be 
shown below. 
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FIG. 1: The present upper and lower limits (95% CL) for 
the Yukawa coupling Xd, for a scalar h (tan/3 if Xv ~ 0, solid 
line) and pseudoscalar A (tan/3, dashed line) as a function 
of mass. Upper limits are from the Yukawa process (ALEPH 
and DELPHI results) and lower limits from the Z h{A)'y 
processes measured at LEP. In addition the upper 90% CL 
constraints from the T decay (K), rescaled by a factor 2 are 
shown. TEVATRON upper limits, both for h and A, are 
presented. The lines denoted "g — 2{B)" corresponds to the 
upper limits for h or A obtained in this paper for case B and 
L'm(95%): in the upper panel one-loop results, in the lower 
one the results of the two-loop analysis (for Xv = 0), see text 
for details. 
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C. 2HDM(II) contribution to 

In the 2IIDM (II) neutral scalars h and H, pseu- 
doscalar A, as well the charged Higgs boson can 
contribute to a^. 

level 



the one-loo 




80 



There exist relevant calculations at 
H, |l| and the two-loop level 
, see also earlier papers [|l5[ In con- 
trast to the one-loop approach where each Higgs boson 
exchange is given by a separate diagram (Fig. 2), var- 
ious Higgs particles may appear in the same two-loop 
diagram, see Fig. 5 (right). 

We assume that the lightest Higgs boson, h or A, dom- 
inates the Ml 2HDM (II) contribution, i.e. aff^^^ « a'^, 
or (a simple approach, see also ||6^ ). This approach 
should hold for masses below 50 GeV, as discussed in Sec. 
III.B. For higher masses, which also are considered here, 
this should be treated as an assumption of a large gap 
between Higgs bosons masses. Since the charged Higgs 
boson mass should be bigger than 320 GeV we do 
not consider to be a lightest particle of the model. 

We calculate separately the one-loop and two-loop con- 
tributions to both and a^. This way the importance of 
the two-loop diagrams can be seen, moreover our results 
can be easily compared with other one-loop calculation , 
e.g. For a pseudoscalar A the two-loop contribu- 

tions are due to diagrams with fermion loops only. For a 
scalar the fermionic, W and a charged Higgs boson loops 
can contribute. However, the W and charged Higgs bo- 
son contributions are expected to be strongly suppressed 
for a small mass of h, where, according to the LEP data, 
the coupling Xv should be small (see discussion in Sec. 
III.B). For a simplicity in calculating of the two-loop con- 
tribution for h we explicitly assume for a whole mass 
range Xv ~ sin(/3 — a)=0, what means a domination of 
the fermionic loops in a''. The more extensive consider- 
ation will be given elsewhere p6| . 

For one- and two-loop approach we derive constraints 
on Yukawa coupling for h and A obtained from the esti- 



mated dUfj, intervals (Sec. II. B), by taking a 



2HDM 



San 



1. One-loop calculation 
The set of the relevant diagrams is presented in Fig. 



with:Qh{x) = x^(2 — x), QAix) = —x^ 
Higgs particle exchange is described by: 



The charged 



'"V ~xn-x) 

z I dx- ^ 



l)z + l 



The scalar contribution a^ |105| is positive whereas the 



pseudoscalar and the charged Higgs boson give negative 
contributions. Each one-loop contribution disappears 
in the limit of large mass like m^/A^f^ ln(M^/m^), see 
also |6^, . At small mass each contribution reaches its 
maximum (or minimum if negative) value. 

The individual (absolute value of) contributions, with 
couplings as in the SM, i.e. with Xd ~ ^ ' shown in 
Fig. 3. For mass above 0.2 GeV the charged Higgs con- 
tribution is much smaller than the contributions due to 
neutral Higgs bosons h and A. One observes that the 
h and A give practically the same contribution for mass 
larger than few GeV, only difference being in sign. In 
principle one should take into account a possible cancel- 
lation of these contributions, especially for masses above 
50 GeV, see ||. 

b. Constraints. First we present constraints based 
on the standard 95% CL intervals, L'm(95%). For case B 
this approach leads to the upper limits on the Yukawa 
coupling Xd for a pseudoscalar (tan /3) and for a scalar 
(tan/3 if /3 = a). Results are shown in Fig. 1 (upper) 
as lines denoted "g — 2(B):one-loop" . For case A one 
obtains an allowed band for the Yukawa coupling of a 
scalar only (only h gives a positive contribution). It can 
be found in Fig. 4 (only lower edge is shown) and in Fig. 6 
(upper), as a region between lines denoted A/B and A. 
A pseudoscalar is excluded in this case. 

Both for B and A case, the limits for Xd rise with mass 
of the Higgs particle, what reflects the decrease of the 
corresponding a'^^'^ terms with increasing mass. For case 
A, only large tan (3 (for P — a) greater than 10 is allowed 
for mass of h above 5 GeV, see Fig. 4. 

Next we discuss results based on lim+{95%) interval. 
For case A one obtains the allowed band for the Yukawa 
coupling Xd^ which practically overlap with the consid- 
ered above lim{95%) band, presented in Figs. 4 and 6(up- 
per). Now, also for case B instead of the upper limits 
an allowed band for the Yukawa coupling appears for h 



a. Individual contributions. The relevant one-loop 
formulae from the Appendix of paper , based on re- 
sults |l|, |l|, |l|l, are given by (A = h, A or iJ±) 

„A| — Jl_T f — ^ "^^ ^A 

Unlone-loop — „ ^ ~ n l\ /r ^d' 

on I M\fi/ 

If /3 = a, the coupling x^ is universal for /i, A and -ff^, 
and it is equal to tan/3 (see sec. HI. A). 

The integral for the neutral Higgs boson contri- 

bution is given by (with z ~ [m^/ Mf^Y): 



Lh(A){z)^z\ dx 



x^z + (1 — x) 




- h,A ' 




IJ'~ 



FIG. 2: One-loop contribution to g — 2 for muon due to a 
neutral scalar h (or H), pseudoscalar A and a charged Higgs 
boson H'^ exchange. 
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Exclusion 95%CL for h in 2HDIVI 



3 

E 
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1 
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0.01 
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contribution to (g-2)/2 for muon 



: 1 ' ' ' 

h;;rJT&A: one-loop 
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A--- ^ 




s^^jMA: two-loop^ 


- H+: ... 




■_ one-loop 




- in10'^{-11} 
1 . . . 


.... 1 1 ■ 



0.1 



1 10 
M_h(A,H+), GeV 



100 



FIG. 3: The (absolute value of) individual contributions 
to from a neutral scalar h (solid line), a pseudoscalar A 
(dashed line) and a charged Higgs boson H'^ (dotted line). 
For results based on one-loop calculation: A and iJ"*" contri- 
butions are negative. Two- loop diagram contributions only 
for A and h (denoted "1") are based on the down-type fermion 
loops. For h also results with additional charged Higgs boson 
loop are shown (coupling of equal to a first term in eq.6. 
with A* = 0) line "2" ("3") corresponds to =800(400) 
GeV. The two- loop h (A) contribution is negative (positive). 
Yukawa couplings as in the SM are assumed. 



(only). In Fig. 6 (upper) these two allowed bands for h 
are compared: a region between lines A/B and A, and a 
wider region between lines A/B and B+. For both cases, 
A and B, one can exclude the negative contribution, as 
they may be realized at the level 0.5 % or 2.8%, respec- 
tively. This means an exclusion of a pseudoscalar. 



C 

CO 





1 1 1 1 1 


1 1 \ \ >/ 


1 uu 


: TUKdWd 

L_ A^ 




10 


\ Upsilon K 


g-2 (lower limit) 


1 


1 





10 

Mh, GeV 

FIG. 4: The exclusion plot for a Yukawa coupling Xd (tan /3 
for fi = a) for a scalar in the 2HDM (II) (a one-loop calcula- 
tion). The lower limit from the (g — 2)^^ data for case A, the 
allowed region lays above the line with arrows. Upper limits 
from the Yukawa process (A=ALEPH and D=DELPHI), and 
the T decay (K). White areas are allowed at 95% CL (both 
lim and lim+). 



h,A 
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h 




H+,W+ 
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FIG. 5: Two-loop contributions to from a light h or A 
with a fermionic loop (left); two-loop diagram for a light h 
with a charged Higgs boson H'^ or loop (right). 



2. Two-loop calculation 



The two-loop diagrams, see Fig. 5, can give large con- 
tributions since they allow to avoid one small Yukawa 
coupling with muon in favor of the coupling with other, 
potentially heavy, particles circulating in the loop |jl^, |l^, 
|80| , |l[ In addition to the mass effect, such contri- 
bution can be enhanced further by an additional factors 
(eqs. 1-3,5). 

In principle diagrams with a Z boson, instead of the 
exchanged photon, may appear. However they are ex- 
pected to be small and are neglected in this analysis. 

Below we consider the diagrams with fermionic loops 
Fig. 5 (left), which are the only two- loop contributions 
for a pseudoscalar A. For h in addition also a charged 



Higgs boson loop (Fig. 5 (right)) is taken into account. In 
calculation of the two-loop contributions we take Xv~^' 
as discussed above. With this condition one can neglect 
a M^-loop contribution. 

a. Individual contributions. The contributions from 
diagrams with fermionic loops presented in Fig. 5 ( left ) 
are given by the following formulae for A — h,A ||8y, |8^, 
flOl: 



'.f,\two-loop - ^K^A' /a - n , Xd- (7) 
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The K parameter is equal to 1 for a pseudoscalar A, for 
a scalar h we have 



h 

K=^ for/ 



' u — type quarks''' , 



(8) 



K = 1 for f = "d — type quarks" 
for / = " charged leptons" . 



The ^ parameter is equal to 1 for leptons and to NcQg 
for a quark q with the electric charge Qqe, Nc = 3. 

The integral -^{(^) for the neutral Higgs boson contri- 
bution with a fermionic loop is given by: 



L 



Q{x)h{A) ^^x{l-x) 



dx 

Q x{l — X) — z 



(9) 



ruf 



Ml 



hi A 



and Qh{x) = —(1 — 2x(l — x)), Qa{x) = 1. A diagram 
with the loop presented in Fig. 5 (right) contributes to 
a^ with the integral Lf^^ (x) given by a similar expression 

as for fermionic loop (eq.ll), with Qf^^ {x) — --(a;(l — a;)) 
prf . The corresponding coefficients for its contribution 
to a'^ are: ^ = 1 and k = XH+/Xd- 

The integrals describing two-loop contribution with a 
fermionic loop is negative for h, while positive for A, on 
contrary to the one-loop results. In Fig. 3 the (absolute 
value of) two-loop fermionic contributions, with /i, r, 5 
loops for h (denoted "1", solid line) and A (dashed line), 
are presented for Xd ~ ^- These two-loop contributions 
dominate for masses above few GeV over the correspond- 
ing (absolute value of) one-loop terms for both h and A. 
They have a milder than the one-loop contribution de- 
pendence on the mass of h (or A), namely ln(m^/M^). 
From the figure one can read that the change of the sign 
of the sum of the one- and two-loop terms appears at 
mass of A (h) equal to 3 (5) GeV. 

The top contribution dominates for Xm = 1 in the large 
mass region, being of order ~ 1-10~^^ (not shown). How- 
ever it does not play an important role in the present 
analysis, since its contribution is proportional to XdXu 
equal to -1 for the pseudoscalar A and (if Xv~^) ^^^'^ 
for h. So, the top-loop contribution has no additional 
enhancement factor, and for large Xd the d-type fermion 
loops dominate, both in ajj and a^. 

The charged Higgs boson loop contribution gives, for 



Xv- 



-0 and K 



MhZlL. 

H + 



> 0, a negative contribution to 



a^. This contribution (absolute value of) rises with a 
mass of h like Ml/ Mjj+{\nMjj+ /Ml + 5/3). Its effect 
can be seen at Mh above 100 GeV. In Fig. 3 these results 
for ajj are presented for = Q (fines "2" for =800 
GeV and "3" (400 GeV)). It is clear that a charged Higgs 
boson loop with the parameters as described above leads 
to a small modification at large mass of h only. More 
detailed discussion will be given elsewhere [Q . 



h. Constraints. A full two-loop calculation, with 
both one- and two-loop diagrams included, leads to re- 
sults which differ significantly from the ones based on the 
one-loop diagrams only. The main difference is related 
to the fact that for masses above 5 (3) GeV, a scalar 
(pseudoscalar) contribution to has opposite sign as 
compared to the corresponding one-loop contribution. It 
means that now for each scenario, with a light h and a 
light A, the contribution can be positive or negative de- 
pending on mass. It is not the case for the one-loop terms 
where the corresponding contributions have a fixed sign. 

In the derivation of the constraints in two-loop ap- 
proach we take also into account the non-leading terms, 
namely the top contribution and for h also a term XuXd 
due to the charged Higgs boson loop (6). The results are 
as follows. 

The case B allows in the standard approach for a neg- 
ative and positive contributions to a^. Therefore the 
upper limits for Xd exist for the whole mass ranges for 
both h and A. In Fig. 1 (lower ) we present the obtained 
limits (lines denoted "g-2(B):two-loop" ) for a mass re- 
gion above 5 GeV. Here the upper limit for pseudoscalar 
was obtained from the maximal positive i5a^ value, while 
for a scalar - from a maximal negative one, see the first 
table in Sec. II. B. This explains why a scalar is now 
constrained more tightly than a pseudoscalar, contrary 
to limits based on the one-loop calculation (Fig. 1 (up- 
per)). For a scalar results obtained with only down- type 
fermions /i, r, b included in the two-loop calculation are 
represented by the line denoted "1" . Lines "2" and "3" 
are obtained if in addition one takes into account the 
charged Higgs boson-loop, as described above. Note that 
if one compares the two-loop constraints for a scalar to 
the one-loop constraints for a pseudoscalar these for a 
scalar are more tight for mass above few GeV, due to a 
weaker mass dependence of the a'/^\two-ioop, see Fig. 3. 

The case A leads to an allowed band for a positive con- 
tribution to a^, i.e. for a scalar h with mass below 5 GeV, 
and for a pseudoscalar A with mass above 3 GeV, sim- 
ilar results were found in |8^, Obtained constraints 
are presented in Fig. 6 (lower) and in Fig. 7 as the re- 
gions between lines A/B and A. In the same figures also 
results for case B obtained with ^im+(95%) are shown 
(regions between lines A/B and B+). It is not clear how- 
ever whether this approach should be used here, as for 
both, a light A and a light h, scenarios, both positive and 
negative contributions to are possible, as discussed in 
Sec. III. C. 



D. Combined 95% CL constraints of the 2HDM 
(11) 

When the above constraints obtained from a new g — 2 
for the muon measurement are added to the existing con- 
straints from other processes discussed in Sec. HI. B. 
(and also in |5^) interesting conclusions can be reached 
in the 2HDM (II) for both scenarios: with a scalar h and 
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a pseudoscalar A being the lightest particles in the model. 
We start discussion of the results based on one-loop cal- 
culation which can be compared with similar analysis 
ifzsl , then results of the two- loop analysis are presented. 

The constraints are obtained for (absolute value of) 
Xdy which is equivalent to tan/3 for pseudoscalar and for 
scalar, provided in the latter case Xv equal or close to 
zero. Only in the calculation of the two- loop contribution 
we explicitly use the assumption Xv= 0- For a simplicity 
of the discussion we will use below the tan /? to represent 
the Yukawa coupling Xd for £^ scalar in all cases. 



1. Allowed regions from one-loop results. 

The upper limits for tan /3 for a pseudoscalar A which 
were obtained from .g — 2 data for case B are much tighter 
than the limits from other experiments for mass above 10 
GeV. Still a window with a light A is open for tan (3 below 
~ 10 (Fig. l(upper)). For a scalar h a weak improvement 
is observed only for a mass range between 60-70 GeV 
where tan/3 has to be lower than ^ 300. 

The {g — 2)f^ results for case A rule out a pseudoscalar, 
for a scalar they improve considerably existing limits. 
The obtained an allowed band in the (tan /?, Mh) plane, 
is equivalent to existance of both the upper (as above) 
and also the lower limits for the Yukawa coupling. An 
allowed by all experiments region appears for a scalar 
h with mass between 40 and 70 GeV at tan/3 above 80 
(Fig. 6 (upper)). In addition a small region of mass of 
h around 10 GeV and tan between 25 and 35 remains 
allowed (Fig. 4), see also |78| where the same result was 
obtained. Note, however, that this small allowed area can 
be closed by taking into account the mentioned in Sec. 
III.B constraint from the global fit ||5^. This constraint, 
for the charged Higgs boson mass above 300 GeV, forb ids 
tan/3 to be larger than 22 for Mh=10 GeV flO^ F08 |. 

At the one-loop level it is sensible to study ^im+(95%) 
intervals, since two scenarios of the 2HDM(II), with a 
light h or A, give contributions of a definite sign inde- 
pendently on mass. Pseudoscalar, giving a negative con- 
tribution, is ruled out. An allowed lim+{95%) band for 
case B obtained for a scalar is wider than for discussed 
above case A (compare the region between lines A/B and 
B+ with the region between A/B and A in Fig. 6 (upper)). 
Therefore also the allowed region in the (tan /3, M^) plane 
is much larger for case B than for case A. It starts at mass 
10 GeV and stops at 70 GeV with tan (3 between 10 and 
~ 300. 



In Fig. 6 (upper) we plot also an the expected upper 95 
% CL limits from process gg h tt at the ep collider 
HERA (dotted lines for lower and higher luminosity), 
from [ |60| , |6l| . These measurements may help to cover the 
low mass region for h and A at the intermediate Yukawa 
coupling tan (3. 



2. Allowed regions from two-loop results. 



The two-loop diagrams give dominant contributions to 
a^^^*^ for mass above 3 GeV (5 GeV) for a pseudoscalar 
(scalar) and these contributions have reversed signs rela- 
tive to the one-loop results. Two- loop analysis for these 
regions (based on a sum of the one- and two-loop dia- 
gram contributions) leads to similar conclusions as the 
one-loop one, with a reversed role of a scalar and a pseu- 
doscalar. In particular now the (standard 95% CL) up- 
per limits obtained for case B for scalar, for Mh > 10 
GeV, are much tighter in comparison with limits from 
other experiments (Fig. 1 (lower)). Still a window with 
a light h is open for tan/5 below 10. For a pseudoscalar 
the improvement is weak , being limited to a mass region 
between 40-70 GeV with tan/3 < 100. 

The constraints obtained in case A in a form of allowed 
bands improve considerably the existing up to now limits 
for both scalar and pseudoscalar since in addition to the 
mentioned above upper limits also the lower limits ap- 
pear (Fig. 6 (lower) and Fig. 7). The allowed by {g — 2)^ 
data bands for a scalar h is situated now below Mh—b 
GeV. This region is however excluded by the T decay 
data (Fig. 6 (lower)). On the other hand a pseudoscalar 
with mass above 10 GeV and tan/3 larger than 20 is still 
in agreement with existing data (see Fig. 7). The TEVA- 
TRON data close practically the region of the mass above 
70 GeV for a case A. 

If for case B the lim^ interval is applied this region of 
mass Ma above 70 GeV is still open. 



3. Discussion 

As (g— 1)^ data favor a positive additional contribution 
it is clear that at one-loop level the most stringent limits 
of the new g — 2 data for the muon can be driven for a 
pseudoscalar, the case A leads even to its exclusion. 

The two-loop calculation leads to a radical change of 
the picture. The two-loop diagrams give the dominant 
contributions for mass above 3 GeV (5 GeV) for a pseu- 
doscalar (scalar) and with a reversed sign as compared 
to the one- loop results. Constraints are now tighter for 
a scalar. 

We point out a role in the present analysis of the low 
energy measurement of the Wilczek process T — > h/A'y 
in closing the window for mass below 10 GeV for both a 
scalar and pseudoscalar. The TEVATRON results close 
a part of the large mass and large tan (3 region. The 
global fit adds an important constraints for h with mass 
around 10 GeV. New measurements of Yukawa process at 
LEP by OPAL and DELPHI which were presented 
recently do not change our qualitatively our conclusions, 
see also below. 



10 



IV. CONCLUSION 

We described the room for a new effects as follows 
from the recent (g — 2)^ measurements and from new 
theoretical estimations of aj^*^. For two SM predictions 
arising from different values of a^°'^, case A (DH) and 
case B (J2000), we evaluate 95% CL intervals for a new 
contribution. They can be used to constrain parameters 
of any model beyond the Standard Model. We show how 
these constraints depend on the size of the a^'"' (case A 
and B), and on the type of limits (lim or lim,+ ). The 
upper limits for the negative and positive contributions, 
obtained if the prediction B and the L'm(95%) are used in 
the analysis, are to be contrasted with the allowed band, 
obtained in case A for a positive contribution only. For 
both A and B cases a negative (5a^ contribution is very 
unlikely: a positive (negative) contribution corresponds 
to 99.5 (0.5) % CL for A, while for case B to 97.2 (2.8) 
% CL. At this level models leading to only negative Sa^ 
can be excluded, and vice versa models which lead to the 
positive contribution only should be accepted as possible 
models. For such models the derived lim+{95%) estima- 
tions differ from those based on lim{95%) considerably 
for case B. 

We applied the obtained intervals to constrain param- 
eters of the 2HDM (II) using a simple approach, where 
only one Higgs boson, h or A, contributes. In the one- 
loop calculation a light scalar scenario leads to the posi- 
tive, whereas the one with a light pseudoscalar to the neg- 
ative contribution to a^, independently of mass. In the 
two-loop analysis, based on a sum of the one- and two- 
loop diagram contributions, the situation changes dras- 
tically. Now the positive contribution can be ascribed 
to a scalar h with mass below 5 GeV or a pseudoscalar 
A with mass above 3 GeV. In both, one- and two-loop, 



approaches we derive tight constraints on the Yukawa 
couplings to muon of h and A. When these constraints 
are combined with constraints arising from other exper- 
iments, especially the Wilczek process and the Yukawa 
processes at LEP and TEVATRON, large part of the pa- 
rameter space for light h or light A can be excluded. 

Our results obtained for the case B lead to an improved 
upper limits for both h and A. For a more constraining 
case A our results including all other existing constraints 
are as follows (in parenthesis the limits obtained if the 
newest DELPHI data are included ): 

The one-loop calculation excludes a pseudoscalar while 
allows for an existance of a light scalar h with mass above 
40 GeV (50 GeV) and below 70 GeV, and tan/3 larger 
than 90. 

The two-loop analysis allows for an existance of a pseu- 
doscalar with mass between ~ 10 GeV (25 GeV) and 70 
GeV, and tan/? above 20 (30). A light scalar being ex- 
cluded by combining the (g — 2)^ and the T decay data. 
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Exclusion 95%CL for h in 2HDM (II) 
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FIG. 6: Upper panel: one-loop result for h. Limits for the Yukawa coupling of the scalar h normalized to the SM value, Xd 
(equal to tan f3 for /? = a) as a function of the mass Mh ■ The constraint from the 77 decay excludes 95% CL the mass below 
0.28 GeV (vertical line), upper 95% CL limits from the Yukawa process (ALEPH measurements for pseudoscalar - dashed line, 
DELPHI data for scalar (pseudoscalar) - solid (dashed) line). Upper 90% CL limits from the T h/A'y, (K, N and L results, 
see text). Results from (K) are rescaled by a factor 2. The upper 95% limits from the TEVATRON collider for h and A 
(rescaled by \/2 for masses below 130 GeV) are also shown (dot-dashed line). Allowed Zjm(95%) bands for h from the newest 
g — 2 data for the muon: regions between the line corresponding to both A and B cases (denoted A/B) and lines giving lower 
bounds calculated for case A and case B (H?n+(95%)) (denoted B+), respectively. For a comparison the expected upper 95 % 
CL limits from process —> /i —> rr at the ep collider HERA are shown (dotted lines for lower and higher luminosity). Lower 
panel: sum of one- and two- loop results for h (for P = a). Curves as in the upper panel. 
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Exclusion 95%CL for A in 2HDM (II) - two loop 
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FIG. 7: Two-loop result for A. Limits for the Yukawa coupling of the pseudoscalar A normalized to the SM value, Xd equal 
to tan/3, as a function of the mass Ma- The allowed bands from g — 2 data for A and B(Zim+ (95%)) cases are shown together 
with constraints from other experiments, details as in Figs. 6. The white area is allowed by the all existing data. 



